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ABSTRACT: Superconducting correlations in the two dimensional t — J model at zero 
temperature are evaluated using numerical techniques. At the fermionic density (n) ~ 
1/2, strong signals of d x 2_ y 2 superconductivity were observed in the ground state. These 
conclusions are based on a study of static pairing correlations, the Meissner effect, flux 
quantization, and other indicators of superconductivity. It is argued that these results can 
be explained using a spin dimer "liquid" state. A phase diagram of the two dimensional 
t — J model is presented. 
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The study of high-T c superconductors continues attracting considerable attention. Re- 
cent calculations suggest that non s-wave symmetry pairing interactions with nodes, may 
explain some of the unusual properties of the cuprate compounds, in particular the be- 
havior of relaxation rates in the YBa2Cu307 material, as well as the systematic presence 
of spectral weight inside the superconducting gap. 1 More specifically, the possibility of 
d x 2_ y 2 superconductivity in the cuprate materials has been recently discussed. 2 ' 3 Most of 
these calculations have been performed without specifying the details of the interaction, 
but analyzing a BCS-like gap equation for different pairing symmetries. Thus, it would 
be important to find a realistic model of strongly interacting electrons having a d x 2_ y 2 
symmetric superconducting state as ground state. From the properties of this state, dy- 
namical responses of a d-wave condensate could be studied, and concrete predictions would 
be made to contrast theory with experiments. In this scenario, numerical studies are im- 
portant to decide whether a given electronic model presents a superconducting phase, 
specially since the strongly interacting character of several realistic models makes most 
analytical approximations questionable. The one band Hubbard model is a typical exam- 
ple of these problems, i.e. while for some time it was assumed that the ground state at 
finite hole doping superconducts, Quantum Monte Carlo simulations have not supported 
these claims. Then, the issue of whether purely electronic models of high-T c materials 
present a superconducting ground state is still open. 

The purpose of this paper is to present numerical results suggesting that the widely 
studied two dimensional t — J model has a superconducting phase in a previously unex- 
plored region of parameter space. The symmetry of the condensate is d x 2_ y 2, and thus this 
model may become a physical realization of the d-wave pairing scenarios recently proposed 
in the literature. 3 ' 2 The superconducting phase observed here appears near the well-known 
region of phase separation of the t — J model. 4 ' 5 The possible presence of superconduct- 
ing correlations near phase separation has been recently discussed in other contexts and 
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theories, ' 7 '° but here the first numerical indications are provided that this phenomenon 
may occur in the ground state of a realistic model of strongly correlated electrons. The 
t — J model is defined by the Hamiltonian, 

H = J £( Si .Sj - ln inj ) - t (cJ, s c j)S + h-c), (1) 
(U) (ij}>s 

where cj g denote hole operators; nj = nj j + n^; and square clusters of N sites with pe- 
riodic boundary conditions are considered. The rest of the notation is standard. Here, 
efforts have been concentrated on the exact diagonalization of 4 x 4 lattices, although 
preliminary results for clusters of 20 and 24 sites are available. It has been repeatedly 
shown in the literature that these cluster sizes are large enough to capture the essential 
qualitative physics of several models of strongly correlated electrons. Besides, no other 
available (and unbiased) numerical technique can handle the involved calculations that 
have been carried out for the t — J model without making assumptions about the prop- 
erties of the ground state. To search for indications of superconductivity, let us define 
the singlet pairing operator Aj = cj^(q +x ^ + cj-x^ ± c i+y,| c i-y,|)' wnere + an d — 
corresponds to extended-s and d x 2_ y 2 waves, respectively, and x, y are unit vectors along 
the axis. The pairing-pairing correlation function C(m) = ^ Xli(^i^i+m)? an d its sus- 
ceptibility Xsiup = J2m C(m) have been calculated (where a = d corresponds to d x 2_ y 2 
wave, and a = s to extended-s wave). () denote expectation values in the ground state, 
which is obtained using the Lanczos method. 

Xsup i s shown in Fig. la as a function of J/t, for several densities. In this study, it 
was observed that the d x 2_ y 2 wave susceptibility dominates, presenting at (n) = 1/2 a 
sharp peak at J/t ~ 3. By analyzing several spin and hole correlations, and following 
other criteria 9 it was verified that the fast decay of x^up after the peak is induced by the 
transition to the phase separated region. Changing the fermionic density, it was observed 
that xtup has its maximum value at (n) = 1/2, as shown in Fig. la. xtup has been also 
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evaluated at (n) = 1/2. This susceptibility peaks at approximately the same position as 
Xsup does, but with a smaller intensity. 10 The pairing-pairing correlations as a function of 
distance are shown explicitly in Fig. lb in the region where the susceptibilities have a sharp 
maximum, i.e J/t = 3.0 and (n) = 1/2. As expected from the behavior of Fig. la, the 
dominant correlation functions at the maximum distance on the 4x4 cluster corresponds 
to d x 2_ y 2 symmetry. On the other hand, the correlations for extended-s operators are 
strong at short distances, but decay rapidly at large distances. For completeness, in Fig. lb 
the correlation corresponding to d xy symmetry is also presented. These correlations seem 
more heavily suppressed than for the d x 2_ y 2 and extended-s channels. In Fig.lc, the d- 
wave pairing correlations are shown at (n) = 1/2, as a function of J/t. Their maximum 
value is obtained at the same coupling where Xsup P ea ks, as expected. 

The symmetry of the ground state (obtained with the Lanczos method) under a ro- 
tation of the lattice in ir/2 has been studied. In the region where superconducting cor- 
relations exists, the ground state is odd under this operation, but it is invariant under 
reflexions with respect to the x and y-axis. Then, the ground state at (n) = 1/2 belongs 
to the Big representation of the C^ v group, usually denoted by d x 2_ y 2, in agreement with 
the previous conclusions studying pairing correlations. However, it is worth noticing that 
the state with the lowest energy in the subspace invariant under rotations and reflexions 
(i.e. s-wave) is close in energy to the ground state. More specifically, at (n) = 1/2 and 
J/t = 3.0, the energy of the d-wave ground state is — 26.413t, the lowest s-wave state has 
energy — 26.127t, while, for comparison, the lowest spin one state in the spectrum carries 
an energy — 25.188t. Then, a scenario where s + id pairing occurs in the bulk limit is not 
excluded, although it is clear that the d x 2_ y 2 correlations seem stronger. 

The pairing correlations found in the two dimensional t — J model suggest the existence 
of a superconducting phase near phase separation. To complete the analysis, it is necessary 
to show that a Meissner effect occurs in that region. Recent progress 12 in the analysis 
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of the superfluid density, D s , using linear response theory allow us to carry out such a 
study using techniques similar to those required to analyze the Drude peak in the optical 
conductivity, a(oj), of strongly interacting electrons. 13 Following Scalapino et al., 12 it can 
be shown that D s is given by 

^^-^E^ l<n|Wq)l ° >|2 ' (2) 

n^O 

where e is the electric charge; the current operator in the x-direction with momentum q 

is given by j x (q) = El,<7 ^(c^+x^ ~ 4+±,a d h^ (~ T ) is the kin etic energy operator 
of Eq.(l); |n) are eigenstates of the t — J Hamiltonian with energy E n (where n = cor- 
responds to the ground state), and the rest of the notation is standard. The momentum 
q = (q x , q y ) of the current operator is selected such that q x = and q y — > 0. The con- 
straint of having an infinitesimal but nonzero q y is necessary to avoid a trivial cancellation 
of D s due to rotational and gauge invariance. 12 On the 4x4 cluster, the minimum value 
of q y is 7r/2, and that is the momentum used in the present analysis. D s given by Eq.(2) 
can be evaluated numerically using the continued fraction expansion technique previously 
used to extract dynamical information from finite clusters. 13 In Fig. 2a, D s is shown as a 
function of J/t for several densities. In good agreement with xiupi the superfluid den- 
sity D s presents a sharp maximum in the neighborhood of phase separation at (n) = 1/2 
giving support to the previous conclusions regarding the existence of superconductivity 
in this model. 14 It is interesting to note that the signal is stronger for lower densities, 
e.g. (n) = 0.25, perhaps due to the higher mobility of pairs in that regime. In the phase 
separated region, D s is small, as expected. 

The resistivity of the model has also been analyzed. From previous studies of the op- 
tical conductivity in strongly interacting models, 13 it can be shown that the Drude peak 
is given by a simple modification of Eq.(2), i.e. it is enough to replace D s — > Drj rU( j e , and 
consider zero momentum, q = (0,0), in the current. 13 ' 12 In Fig. 2b, D Dru( j e is shown as 
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a function of J/t, for several densities. In the region of phase separation, the conductiv- 
ity is small as expected, while for smaller values of J/t, the Drude peak is considerably 
larger. A finite value of Drj ruc j e m the bulk limit implies a zero resistivity, p = 0, since 
a(u — > 0) = a? -1 = DDrude^t^)- Fig. 2b suggests that this result will hold not only in the 
superconducting region, but it will survive a further reduction of the coupling into the 
small J/t regime, i.e. even in a phase without pairing. This example shows that p is not 
enough to distinguish between a "perfect metal" and a "superconductor", and thus the 
previously discussed study of the superfluid density is crucial to show unambiguously the 
presence of superconducting correlations in the model. 12 To further complete the present 
analysis, the response of the system to an external magnetic flux <p was studied. For this 
purpose, a phase factor e^l N is introduced in the kinetic energy hopping terms of Eq.(l), 
but only in the x-direction. This is equivalent to allowing a nonzero flux across one of the 
"holes" of the torus. 15 In Fig. 3a, the ground state energy AE(<p) = E(<p) — E(<p = 0), in 
the zero momentum subspace, is shown as a function of 0, at density (n) = 1/2. In the 
region of pairing, J/t = 3.0, the energy presents two minima, one located at cj) = (mod 
2tt), and a nontrivial one at = tt, signaling the presence of carriers with charge 2e in the 
ground state, in agreement with the analysis based on the pairing correlations. 

What is the nature of the superconducting state at (n) ~ 1/2? It is reasonable to expect 
that the same force that produces phase separation, is responsible for superconductivity. 
Actually, if two electrons are considered on an otherwise empty lattice, they form a bound 
state at J/t = 2, and at low electronic density this same attraction leads to phase separation 
when the coupling is increased. 4 In this respect, the antiferromagnetic coupling should be 
considered as an attractive interaction in the Hamiltonian at low densities, and thus the 
presence of superconductivity in the model is easily understood. At large J/t, pairs of 
electrons (at least for small (n)) are expected to have a size comparable to the range of the 
force, namely approximately one lattice spacing, and in this respect the superconducting 
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correlations discussed in this paper may be of the bipolaronic type. The pairs should be 
coupled in short spin singlets forming dinners. To check these ideas let us analyze the spin- 
spin correlations. In this scenario, each electron is coupled with only one other particle in a 
spin singlet, but due to rotational invariance, that particle can be located at any of the four 
possible nearest neighbors. Then, the correlation at distance of one lattice spacing, should 
be 1/4 of the on-site correlation, and it should vanish at larger distances. The results 
shown in Fig. 3b obtained at (n) = 1/2 and J/t = 3.0, are in excellent agreement with 
this picture. Another issue to address is the possible formation of a "crystal" structure. 
Is there any special order in the position of these dimers? For that purpose hole-hole 
correlations, h(m) = (nj 1 (0)n^(m)) (where n^ is the hole number operator) were studied, 
i.e. once a hole is located at a given site 0, then correlations with other holes are evaluated. 
Asymptotically, h(m) should decay to (n) at large distance. In Fig. 3c, h(m) is shown at 
(n) = 1/2 and J/t = 3.0. The hole-hole correlations rapidly decay to its asymptotic 
value, showing that there is no special pattern in the hole distribution (or, equivalently 
at (n) = 1/2, in the electronic distribution). The spin-gap in the neighborhood of phase 
separation has also been studied. Although to obtain conclusive results a careful finite size 
scaling analysis is necessary, the data suggest the presence of a finite spin-gap, 16 compatible 
with the idea of small size dimers as responsible for the superconducting correlations. 17 
To summarize these ideas, in Fig. 4a a snapshot of the "dimer liquid" state that here is 
claimed to be compatible with the numerical results is presented. 18 Note that this state 
explains the observed maximum that Xsup presents at (n) = 1/2, since at that density the 
model has the maximum number possible of mobile dimers. At smaller densities, there are 
fewer dimers contributing to the signal; while closer to (n) = 1, the pairs of electrons have 
less mobility for lack of space. A similar state, although defined in a rigid crystal pattern, 
was previously discussed in the context of the t — J — V model. 7 ' 6 

Summarizing, in this paper numerical evidence suggesting that the t — J model in two 
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dimensions has a superconducting phase at zero temperature has been discussed. The 
pairing correlations are the strongest at density (n) = 1/2, and near phase separation. 4 ' 5 
The symmetry of the pairing state corresponds to d x 2_ y 2. The Meissner effect, as well as 
flux quantization calculations support this scenario. The size of the pairs seem small at 
(n) = 1/2 as suggested by the spin-spin correlations, and they are in a disordered state. 
Thus, a liquid of dimers may represent the physics of this condensate. These results have 
several implications: i) they are the first numerical evidence that the t — J model super- 
conducts in two dimensions. Previous numerical studies 13 concentrated their efforts near 
(n) ~ 1, but in that regime the signal for superconductivity would be too weak to be 
detectable; ii) In addition, it was found that the symmetry of the superconducting con- 
densate is d x 2_ y 2, and thus this model may become a realization of recent proposals to 
explain the phenomenology of high-T c materials making use of non s-wave pairing interac- 
tions with nodes. 2 ' 3 Based on the present calculation and others 4 ' 5 the currently available 
information for the phase diagram of the two dimensional t — J model at zero temperature 
is sketched in Fig. 4b. The notation is explained in the caption. The "binding" region 
denotes a regime where pairs are formed, but they are not condensed in a superconducting 
state. The details of this phase diagram close to half-filling are more difficult to address 
numerically than at (n) = 1/2. However, the possibility that the model superconducts 
also at low hole doping is not excluded. Whether there is an analytical continuation be- 
tween (n) = 1/2 and large J/t, and densities closer to half-filling and smaller couplings is 
a crucial issue for the success of the t — J model as a phenomenological model of high-T c 
superconductors. This important subject will be addressed in future publications. 

This work benefited from useful conversations with A. Moreo, M. Luchini, D. Scalapino, 
R. Laughlin, and F. Ortolani. The authors thank the Supercomputer Computations Re- 
search Institute (SCRI) and NCSA, Urbana, Illinois, for their support. 
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Figure Captions 



la d x 2„ y 2 superconducting susceptibility, Xs up , as a function of J/t, at densities (n) = 0.25 
(A), (n) = 1/2 (□), and (n) = 0.75 (□). 

lb Pairing-pairing correlation function C(m) as a function of distance m, at density (n) = 
1/2 and J/t = 3.0. □ denotes d x 2_ y 2 pairing correlations, A indicates extended s 
correlations, while □ corresponds to d xy correlations. 

lc Pairing-pairing correlation function C(m) in the d x 2_ y 2 channel, as a function of dis- 
tance m, at density (n) = 1/2. A, □ and □ are results for J/t = 1.0, 3.0 and 4.0, 
respectively. 

2a Superfluid density, D s , versus J/t, at several fermionic densities. □ corresponds to 
(n) = 1/2, A denotes results for (n) = 0.25, while □ indicates (n) = 0.75. 

2b Drude peak, Dj3 ru( j e , as a function of J/t for several densities. The notation is as in 
Fig.2a. 

3a Energy of the ground state as a function of an external magnetic flux (p. The energy 
at zero flux is subtracted from the result i.e. AE((/)) = E((/)) — E(0). The subspace 
of zero momentum is considered, and the density is (n) = 1/2. □ denotes results at 
J/t = 3.0, while □ corresponds to J/t = 4.0 i.e. inside the phase separated region. 

3b Spin-spin correlation S(m) as a function of distance at density (n) = 1/2, and coupling 
J/t = 3.0 i.e. in the superconducting region. 

3c Hole-hole correlations h(m) as a function of distance at density (n) = 1/2, and coupling 
J/t = 3.0 i.e. in the superconducting region. 

4a Qualitative representation of a "dimer liquid" state, presumed to be the ground state 
in the superconducting region discussed in this paper. 
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4b Schematic semi — quantitative phase diagram of the t — J model in two dimensions 
at zero temperature, as a function of coupling J/t, and hole density x = 1 — (n). The 
curves separating the region at small J/t, presumably a Fermi liquid, FL, from the 
"binding" region, as well as the separation between binding and d x 2_ y 2 wave super- 
conductivity, are rough estimations based on the study of binding energies, and the 
strength of Xsup- The transition leading to phase separation is more accurate, and 
in qualitative agreement with high temperature expansions. 5 Near half-filling, the cal- 
culations are more difficult, and it is only known that antiferromagnetic, AF, and 
ferromagnetic, FM, correlations are important. 
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